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We investigate flow properties of immiscible Bose-Einstein condensates composed of two differ¬ 
ent Zeeman spin states of ®^Rb. Spatially overlapping two condensates in the optical trap are 
prepared by application of a resonant radio frequency pulse, and then the magnetic field gradient 
is applied in order to produce the atomic flow. We find that the spontaneous multiple domain 
formation arising from the immiscible nature drastically changes the fluidity. The homogeneously 
overlapping condensates readily separate under the magnetic field gradient, and they form stable 
configuration composed of the two layers. In contrast, the relative flow between two condensates 
are largely suppressed in the case where the magnetic field gradient is applied after spontaneous 
domain formation. 

PACS numbers: 67.85.-d, 03.75.Kk, 03.75.Mn, 05.30.Jp 


I. INTRODUCTION 

Eormation of heterogeneous spatial structures in bi¬ 
nary phase-separating fluids and their fluidity are in¬ 
triguing objects of research in the various fields such as 
non-equilibrium physics, fluid mechanics, and rheology 
[l|. Multi-component atomic Bose-Einstein condensates 
(BECs) realized in recent decades are quantum fluids 
that display phase separation Q. The unique feature 
of atomic BECs is the unprecedented controllability of 
relevant experimental parameters, which provides novel 
approaches for studying the physics of phase separation. 
Eor example, atomic interaction strengths can be tuned 
via Eeshbach resonances Si, and trapping potentials 
can be designed by using optical and magnetic fields. To 
date, miscibility and immiscibility of binary BECs were 
controlled by changing the atomic interaction strengths 
Si and by using Rabi coupling [3, [![. Additionally, 
quantum tunneling across immiscible BECs separated 
into two layers was observed by application of a state de¬ 
pendent potential [9[. Theoretically pattern formations 
at fluid interfaces have been predicted in two-component 
BECs, which are due to interface instabilities such as 
Rayleigh- Taylor [isini and Kelvin-Helmholtz instabili¬ 
ties [13, [ 3 . 

When two spatially overlapping immiscible BECs are 
prepared, domain structures are spontaneously formed 
so that the free energy of the whole system is minimized 
[iM^. This spontaneous spatial structure formation 
has been widely observed in dual species BECs [5[ and 
two-component BECs composed of two different hyper- 
fine states m, [13, Ell or Zeeman states [H, Ell, EqI . In this 
paper, we investigate influences of spontaneous structure 
formation on the fluidity of immiscible two-component 
BECs, and conversely influences of non-equilibrium fluid 
flow on the domain structures. Spatially overlapping two- 
component BECs comprised of different Zeeman states of 
^^Rb atoms are prepared by application of a resonant ra¬ 
dio frequency (rf) pulse. The two-component BECs ho¬ 


mogeneously overlap just after the rf pulse, and domain 
structures grow distinctly after 40 ms. In order to gener¬ 
ate the relative flow between the two-component BECs, 
a spin-dependent potential gradient is applied. We find 
that the center of mass (c.m.) motion is drastically 
changed depending on whether the two-component BECs 
are homogeneously overlapped or multiple domain struc¬ 
tures are formed. When two-component BECs are homo¬ 
geneously overlapped, they readily separate into two lay¬ 
ers, which is a stable configuration for immiscible fluids. 
In contrast, multiple domains refuse to separate against 
the spin-dependent potential gradient, and the relative 
flow is suppressed. 

In the previous experiments in Refs. [23, [111? 
collective dynamics of immiscible BECs in a highly 
displaced magnetic harmonic trap have been observed. 
The effects of periodic collisions between two-component 
BECs on c.m. oscillations and shape oscillations were in¬ 
vestigated. On the other hand, our objective is to reveal 
the influence of spontaneous domain formation on the 
collective motion by controlling the timing of application 
of a spin-dependent gradient potential. 

The paper is organized as follows. In Sec. II, we de¬ 
scribe our experimental and theoretical methods. The 
results and discussions are provided in Sec. Ill, and con¬ 
clusions are given in Sec. IV. 


II. EXPERIMENTAL AND THEORETICAL 
PROCEDURES 

A BEC of 4 X 10^ ^^Rb atoms in the hyperfine 
state |F = 2,mi? = —2) is produced in a crossed far 
off-resonant optical dipole trap (EORT) with axial {z- 
direction) and radial frequencies of (27r) = 32 Hz and 
Cc;r/(27r) = 135 Hz. A more detailed description of cre¬ 
ation of the BEC is given in Ref. [23| . In order to create 
the stable bias field along the z direction, the whole 
apparatus to create the BECs is installed inside a mag- 
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netic shield room, and a laser diode source with low ripple 
noise of less than 20 /iA (Newport, LDX-3232-100V) is 
used as the current source for the z-axis Helmholtz coil. 
This ripple noise current is 50 dB less than the current 
used to create Bz. 

Two-component BECs composed of 

|F = 2 , mi? = — 1 ) = | 1 ) and |F = 2 ,mi? = — 2 ) = | 2 ) 
are prepared by application of a resonant rf 
pulse. The intra- and inter species 5 -wave scat¬ 
tering lengths of the components | 1 ) and | 2 ) are 
(ail, a 22 , <^ 12 ) = (95.68, 98.98, 98.98)aB in units of the 
Bohr radius [26| , which satisfy the immiscible condition, 
y^aiia 22 < CLi 2 ‘ After a holding time of Thoid = 0 ms or 
40 ms, the magnetic field gradient, dBzjdz^ is applied, 
which creates the mi?-dependent potential gradient 
along the 2 ; direction. The two-component BECs move 
under the mi?-dependent gradient potential during Tgrad 
and then they are released from the EORT. The two 
component BECs are then spatially separated along the 
2 ; direction by the Stern-Gerlach (SG) method. After a 
time-of-flight (TOE) of 15 ms, the atomic distribution of 
each component is measured using absorption imaging. 

In the numerical simulations, we use the mean-field 
approximation. The dynamics of the macroscopic wave 
functions Ti(r,t) and T 2 (r,t), which correspond to 
the components | 1 ) and | 2 ), are described by the two- 
component Gross-Pitaevskii (GP) equation: 










( 1 ) 


where m is the mass of a ^^Rb atom, and gj^k = Qkj is 
the coupling constant between \i) and \k). The harmonic 
potential and the magnetic field gradient are represented 
hy Uj{r,t) = m[oj'^{x'^+y‘^)+ulz‘^]/2-jidB{dBz/dz)d{t- 
^hoid)^/2, where ^(•) is the Heaviside step function and 
/iB is the Bohr magneton. In order to include the effects 
of the spin-dependent two-body inelastic loss for mi? = 
— 1 component, we use gi^i = 47 rfi^aii/m — 3 ^/ 162 /14 and 
gi 2 = 922 = 47 rfi^a 22 /^, where 62 = 24.3 x lQ~^^cmVs 
is the inelastic collision coefficient reported in Ref. |^ . 

We solve the two-cornponent GP equation using the 
pseudospectral method The initial conditions are 

prepared as follows: we first obtain the ground state T(r) 
of the GP equation ([ 1 ]) under the condition Ti(r,t) = 0 
by the imaginary-time evolution method, where T(r) 
is normalized by the number of atoms, /dr|T(r)p = 
4 X 10^. The real-time GP equation is then solved for 
the initial condition Ti(r,t = 0 ) = ^2/3T(r) and 
= 0) = ^l/3T(r). In the TOE expansion, the 
harmonic potential in Uj{r,t) is switched off. The nu¬ 
merical meshes are typically 64 x 64 x 512 for in-trap 
simulations and 512 x 512 x 1024 for the TOE simula¬ 
tions. 



FIG. 1: (Color online) rf spectroscopy. The initial state |2) 
is irradiated by an rf pulse. The rf pulse shape is a Gaussian 
function with a standard deviation of 9.5 gs. (a) and (b) show 
the frequency and intensity dependence of population of each 
tjif component, respectively. In (b), the frequency of rf pulse 
is set to be 8.145 MHz. Each point represents the average over 
three measurements with the error bars giving the standard 
deviation over those measurements. 


III. EXPERIMENTAL AND THEORETICAL 
RESULTS 

A. rf spectroscopy and state preparation 

In order to experimentally prepare two-component 
BECs comprised of | 1 ) and | 2 ), we carried out rf spec¬ 
troscopy. Eigure 1 shows the population of each mi? com¬ 
ponent obtained by the application of a single rf pulse 
with a constant frequency to the initial state |2). Erom 
the results of frequency dependence [Eig. 1(a)], the reso¬ 
nant frequency of the rf transition between | 1 ) and | 2 ) is 
found to be 8.145 MHz. The two-photon transition into 
the mi? = 0 component was also observed around 8.133 
MHz. Atoms in mi? = 0, -hi and 4-2 components are 
not populated for the resonance frequency between | 1 ) 
and 12) due to the large second-order Zeeman shift and 
the stable magnetic field environment. As shown in Eig. 
1 (b), we can prepare the superposition state of | 1 ) and 
12 ) with an arbitrary population by changing the inten¬ 
sity of rf pulse at 8.145 MHz. Erom these spectroscopic 
measurements, the strength of the Bz held is found to be 
Bz = 11.685 G. 


B. Spontaneous multiple domain formation in an 
immiscible two-component BECs 

We observed the spontaneous spatial structure forma¬ 
tion of two-component BEGs composed of | 1 ) and | 2 ). 











3 


(a) Experiment 

^hoid = 0 ms 




I I I I I I I 


Z 


| 1 > 

X 


4 



20 ms 40 ms 



60 0 -60 60 0 -60 


60 ms 



80 ms 



. . 



I I I I I I I I I I I I I I 


60 0 -60 60 0 -60 (urn) 


(b) Simulation 



60 ' 'o' -60 60 ' 'o' -60 60 ' 'o' -60 60 ' ' 0 ' -60 60 ' ' 0 ' -60 (^m) 



0 20 40 60 80 

Time after the preparation of two-component BECs (ms) 


FIG. 2: (Color online) Experimentally observed (a) and nu¬ 
merically simulated (b) phase-separation dynamics. Absorp¬ 
tion images are obtained at Thoid = 0 ms, 20 ms, 40 ms, 60 
ms, and 80 ms from left to right. The two-component BECs 
of |1) and 12) states are initially prepared in the ratio of 2:1. 
The magnetic field gradient to create the state-dependent po¬ 
tential is not applied in this experiment and numerical simu¬ 
lation. 


Figure 2(a) shows the typical absorption images observed 
for Thoid = 0 Ills, 20 ms, 40 ms, 60 ms, and 80 ms, where 
the initial ratio of the components |1) and |2) is 2:1. 
Just after the application of the rf pulse (Thoid = 0 ms), 
the two-component BECs overlap each other [leftmost in 
Fig. 2(a)]. After Thoid = 40 ms, spatial domain structure 
grows distinctly [middle and right in Fig. 2(a)]. 

The numerical simulation of the coupled GP equation 
is shown in Fig. 2(b). Though the growth of multi¬ 
ple domain structure is reasonably reproduced in the nu¬ 
merical simulation, the experimentally observed patterns 
are asymmetric along the z-axis. A possible cause is the 
asymmetry of the trapping potential due to experimen¬ 
tal imperfections such as small optical misalignments. In 
addition, the experimentally observed spatial structures 
are not clear due to the finite spatial resolution of our 
imaging system ( 7.5 i^m). 


C. Center of mass motion under the state 
dependent potential gradient 

Figure 3 shows the experimental and numerical re¬ 
sults of the relative c.m. motion between two-component 
BECs, Z 2 — zi, with and without application of the mag¬ 
netic field gradient. In order to investigate the effect 
of the rui^-dependent potential gradient oc —mpdBzIdz 
on two-component BECs without the domain structure, 
we measure the dynamics in two different cases: (i) 
the free evolution without magnetic field gradient (the 


EIG. 3: (Color online) Relative c.m. motion between the 
two-component BECs, — z\^ as a function of time after 
the preparation of two-component BECs. The open circles 
represent the data calculated from results of Eig. 2, namely 
they correspond to the free evolution without magnetic field 
gradient. The filled circles (cross marks) are measured by 
changing Tgrad with the fixed value of Thoid = 0 ms (Thoid = 40 
ms), where the magnetic field gradient of 15 mC/cm is applied 
during Tgrad- The c.m. positions after TOE differ from those 
in the trap, and their values might be affected by the in-trap 
velocity. However the red points in (b) seem to scale those of 
in-trap [Inset in (b)] thanks to the short TOE time of 15 ms 
and small magnetic field gradient of 15 mC/cm. 


open circles) and (ii) application of magnetic field gra¬ 
dient of 15 mC/cm just after the preparation of two- 
component BECs (the filled circles). We define that 
Tgrad and Thoid are the time periods with and without 
the magnetic field gradient, respectively. Thus, Thoid is 
varied with Tgrad = 0 ms in (i), and Tgrad is varied with 
Thoid = 0 ms in (ii). In both of experiments and sim¬ 
ulations, the relative c.m. of the two-component BECs 
is almost unchanged for the case (i), whereas the c.m. 
of the two-component BECs is clearly separated for the 
case (ii). In (ii), although the field gradient pushes the 
two-component BECs in the same direction, component 
|1) moves in the direction opposite to the gradient force. 
This is because component |1) is pushed away by compo¬ 
nent 12) due to the repulsive interaction between the two 
components [^, [^. After Thoid — 40 ms, the relative 
c.m. decreases due to the harmonic potential. This pe¬ 
riod is longer than the free oscillation period in the trap 
due to the inter-component interaction. 

In order to investigate the effect of the domain struc¬ 
ture on the c.m. motion, we observed the dynamics for 
the case (iii): after Thoid = 40 ms the magnetic field 
gradient of 15 mC/cm is applied during Tgrad (the cross 
marks in Fig. 3 are obtained by varying Tgrad)- Since the 
spatial domain structure is observed clearly at Thoid = 
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FIG. 4: (Color online) Structural change and c.m. motion induced by the application of magnetic field gradient, where the 
distance traveled by SG measurement is compensated, (a) Typical absorption images obtained at Thoid = 0 ms and Tgrad = 10 
ms. (b) Typical absorption images obtained at Thoid = 40 ms and Tgrad = 10 ms. (c) c.m. positions (zi and Z 2 ) calculated 
from data of (a) and (b) versus the applied magnetic field gradient. Top and bottom panels represent the experimental and 
numerical results, respectively. In the experimental data, the magnitude of the effective magnetic field gradient is larger than 
the horizontal axis in (c) due to the overshoot of the applied magnetic field gradient. See Supplemental Material for movies 
showing time evolution of the two-component BECs in a trap potential 


40 ms as shown in Fig. 2(a), the magnetic field gradi¬ 
ent is applied after the domain structure is established 
in the case (hi). As well as the case (i), the relative c.m. 
positions are almost unchanged for the case (hi). This 
result indicates that the spontaneous domain formation 
suppresses the relative motion between two-component 
BECs. The quantitative differences between experiments 
and simulations in Fig. 3 are likely due to the asymmet¬ 
ric pattern formation observed in Fig. 2, because the 
quantitative differences largely occur after 40 ms. 

Figures 4(a) and 4(b) show the typical absorption im¬ 
ages obtained for various magnitudes of applied magnetic 
field gradient. When Thoid = 0 ms and Tgrad = 10 ms 
[Fig. 4(a)], in which the magnetic field gradient is applied 
just after the preparation of the two-component BECs, 
two-component BECs are readily separated. The result¬ 
ing spatial structures consist of two layers, which is a sta¬ 
ble configuration for the immiscible fluids under a spin 
dependent gradient potential. In the case of Thoid = 40 
ms and Tgrad = 10 [Fig- 4(b)], even for a large gradi¬ 
ent about 500 mG/cm, the two-component BECs move 
while maintaining multiple domain structure so that the 
spatial overlap between the components |1) and |2) is not 
increased. 

In order to quantitatively evaluate these motions, 2^1 
and Z 2 are calculated from Figs. 4(a) and 4(b), which 


are shown in Fig. 4(c). Compared with case of Thoid = 0 
ms (the open and filled circles), zi moves faster and Z 2 
moves slower in the case of Thoid = 40 ms (the open 
and filled triangles). This result indicates that the mul¬ 
tiple domain structure suppresses the separation of the 
two-component BECs, even when the magnetic field gra¬ 
dient is 500 mG/cm. In other words, one of the two- 
component BECs serves as the potential barrier for the 
other BEC in the presence of the multiple domains. How¬ 
ever, the domain structure does not completely suppress 
the change in 2^1 — 2 ^ 2 . According to the simulations , 
for our experimental parameters, this is because the two- 
component BECs pass by each other, rather than the 
quantum tunneling [9|. In the case that the magnetic 
field gradient of 500 mG/cm is applied from Thoid = 40 
ms, after Tgrad ^ 00 ms, two-component BECs are spa¬ 
tially separated into two layers and the domain structure 
disappears. 

Finally we discuss the fluidity under the larger mag¬ 
netic field gradient. When the typical interaction energy 
^ 12(^1 +ri 2 ) and Zeeman energy fiB{dBzldz)wdomainl‘^ 
are comparable, the effects of the domain become less 
important, where rii and 77.2 are the particle densities of 
the components |1) and |2), respectively, and rcdomain is 
a typical domain size 10 ^ 20 /rm. According to this es¬ 
timation, we find that the effects of the domain becomes 
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irrelevant when dBzjdz ^ 2000niG/cni. 

IV. CONCLUSIONS 

We have investigated the flow properties of immiscible 
two-component BECs composed of two different Zeeman 
states of ^^Rb. The two immiscible BECs evolve in the 
optical trap forming multiple domain structures. In order 
to investigate the effect of these multiple domains on the 
fluidity, a spin dependent gradient potential, which pro¬ 
duces an atomic flow, was applied at different times. We 
find that the induced flow is significantly affected by the 
spatial structures. At early times, the spatially overlap¬ 
ping condensates readily separate under the state depen¬ 
dent potential gradient. However, after sufficient time 
is given to allow domain formation, the two-component 
BECs are difficult to separate to each other, since the 
multiple domains play role of barriers that prevent the 


BECs to separate. Our results demonstrate the com¬ 
plex nature of flow in two-component immiscible BECs. 
Despite the intuitive expectation that immiscible com¬ 
ponents should separate under a spin dependent gradi¬ 
ent potential, multiple domains arising from the immis¬ 
cible nature actually resist the relative flow of the two- 
component BECs. We anticipate that the method we 
have used here to probe the effects of domain formation 
can be applied to further studies of quantum fluids in the 
future. 
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